Phospholipid demixing and the birth of a hpid droplet 



J. Zanghellini0 F. Wodlci, and H. H. von Griinbcrg 
Institute of Chemistry, University of Graz, Heinrichstrafle 28, A-8010 Graz, Austria, EU 

The biogenesis of lipid droplets (LD) in the yeast Saccharomyces cerevisiae was theoretically in- 
vestigated on basis of a biophysical model. In accordance with the prevailing model of LD formation, 
we assumed that neutral lipids oil-out between the membrane leaflets of the endoplasmic reticulum 
(ER), resulting in LD that bud-off when a critical size is reached. 

Mathematically, LD were modeled as spherical protuberances in an otherwise planar ER mem- 
brane. We estimated the local phospholipid composition, and calculated the change in elastic free 
energy of the membrane caused by nascent LD. Based on this model calculation, we found a gradual 
demixing of lipids in the membrane leaflet that goes along with an increase in surface curvature at 
the site of LD formation. During demixing, the phospholipid monolayer was able to gain energy 
during LD growth, which suggested that the formation of curved interfaces was supported by or 
even driven by lipid demixing. In addition, we show that demixing is thermodynamically necessary 
as LD cannot bud-off otherwise. 

In the case of Saccharomyces cerevisiae our model predicts a LD bud-off diameter of about 13 
nm. This diameter is far below the experimentally determined size of typical yeast LD. Thus, we 
concluded that if the standard model of LD formation is valid, LD biogenesis is a two step process. 
Small LD are produced from the ER, which subsequently ripe within the cytosol through a series 
of fusions. 

Keywords: lipid droplet; biogenesis; demixing; organelle dynamic; 



I. INTRODUCTION 

Lipid droplets (LD) arc depots for neutral lipids (NL). 
They exist in virtually all kind of living cells, from bacte- 
ria, to yeasts, to plants and mammals. A LD consists of a 
hydrophobic, NL-containing core surrounded by a phos- 
pholipid (PL) monolayer containing a small amount of 
proteins [l| . In Saccharomyces cerevisiae, LD are mainly 
formed by triacylglycerols (TAG) and steryl esters (SE) 
in roughly equal amounts 0, [1, Conversion of free 
fatty acids (FA) and sterols to NL and their subsequent 
storage in LD is an organism's strategy to risklessly save 
intrinsically toxic FA and sterols for later use. If required, 
FA and sterols may be rapidly released from LD and used 
as pre-fabricated building blocks for membrane lipid syn- 
thesis as well as other complex lipids, and/or as source of 
chemical energy Also, LD are assumed to have a 

function in transporting sterols to the plasma membrane 
Q. Indeed, it is now recognized that rather than being 
inert storage pools, LD are remarkably flexible, dynamic 
organelles jl, [1, Q . 

LD biogenesis is everything but clear. According to a 
widely accepted model, NL accumulate within distinct re- 
gions of the endoplasmic reticulum (ER) membrane, ini- 
tially forming a lens-shaped and then a spherical bulge in 
the membrane (Fig.[T]). After reaching a critical size, ma- 
ture LD will bud-off, being encapsulated in a PL mono- 
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FIG. 1: LD formation according to the standard model. Dur- 
ing their biosynthesis NL "oil-out" in between the leaflets of 
the ER bilayer, forming spherical structures. Mature LD then 
bud-off and form independent organelles. 



layer that is directly derived from the cytoplasmic ER 
leaflet [H, i, i, i, [13 ■ This budding model is in line with 
several experimental findings, (i) In yeast the same ER 
proteins, except for very few, are also detected on the 
LD membrane [ll], • (ii) Most LD proteins lack trans- 
membrane spanning domains p^ . (iii) In yeast mutants, 
which are unable to synthesize TAG and SE, LD do not 
form. Nonetheless all typical LD proteins are found in 
these strains, but now solely localized to the ER [3] and 
the cytosol. It has also been hypothesized that devel- 
oping LD might not bud-off, but are cut out from the 
ER in form of bicelles, leaving a transient hole in the ER 
membrane [isj . 

A nascent LD trapped within the leaflets of the ER has 
never been observed experimentally, thus other mecha- 
nisms for its formation have been suggested. (See [l6l.[T7| 
for reviews.) The most prominent alternative is based on 
vesicular budding [ll, [l^, . In such a process small 
bilayer vesicles are formed, which are subsequently filled 
with NL pn . [2^ . As yet, neither model has been con- 
clusively verified experimentally. Whether this has to be 
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attributed to low resolution of microscopic approaches or 
indicates that the proposed scenarios are wrong is still a 
matter of debate. In this article we take up the former 
position. Moreover, we here assume that LD formation 
takes place according to Fig. [TJ That is to say, NL-fiUed 
bulges are formed in the cytosolic monolayer of the ER 
membrane, from where they subsequently bud-off. Tak- 
ing this scenario seriously, we calculate biophysical con- 
sequences for the process of LD biogenesis. 

Our analysis has been motivated by experimental find- 
ings in yeast that the PL monolayer cornposition of LD 
differs from the one of the ER membrane 0, [H, [l^l . This 
suggests that the PL composition is related to the local 
curvature of the monolayer in agreement with similar ob- 
servation from theory and experiment [25l . [26l |. Lipids 
with a cone-like molecular shape (inducing a positive or 
convex curvature) are expected to be more adapted to 
the spherical surface of LD than wedge-like shaped ones, 
which induce a negative or concave curvature. Thus, the 
former should be enriched on LD surfaces (relative to 
their value in the ER membrane) while the latter ought 
to be depleted. For instance, due to their geometrical 
shape we expect to find more lysophosphatidylcholinc 
(LPC) and less phosphatidylethanolaminc (PE) in the 
membrane of LD compared to the ER. 

In contrast to the well established view that PL demix- 
ing is curvature-dependent, we here argue the converse, 
i.e., that the demixing of a lipid membrane supports the 
generation of membrane curvature. We use yeast LD as 
an example and present calculations based on the stan- 
dard model of LD formation. They reveal that due to 
lipid demixing, the PL monolayer is able to gain energy 
during its shape transition, thus supporting LD forma- 
tion. Additionally, we show that depending on its volume 
and the curvature of its surface, a nascent LD is pre- 
vented from budding-off through a demixing-controUed 
energy barrier, whose height decreases with increasing 
LD volume. We predict that at a LD diameter of about 
13 nm, this barrier completely vanishes and the LD is re- 
leased from the ER. We thus address two essential ques- 
tions 0] : How do nascent LD bud-off from membranes? 
How is the mature size of LD determined? 



II. THEORETICAL MODEL 

To study local deformations in an extended lipid mono- 
layer we use the Helfrich Hamiltonian [l^] , 

C/ = ^dA|^[i7(r)-2co]Vi;A-(r)|, (1) 

which relates the local, total and Gaussian curvature, 
H = ci + C2, and K = ciC2, respectively, of a 2- 
dimensional surface. A, to its elastic energy, U. Here, ci 
and C2 denote the two principal curvatures at any given 
point r on the surface. The thin lipid membrane is char- 
acterized by its spontaneous curvature, cq, and its mean 
and Gaussian bending modulus, k^, and kc, respectively. 



Eq. ([T]) describes a continuous membrane, without any 
reference to its internal structure. However, to be able to 
model a mixed membrane (i.e. a membrane consisting of 
various different types of PL) we assume that ^ is valid 
not only for a monolayer as a whole but also for every 
single PL molecule. Thus, 

U^{r) = Ap,, 1^ [H{r) ~ 2co,f +Kj<{r)'^ , (2) 

denotes the potential energy of a single lipid where the 
index i in the following distinghuishes between the dif- 
ferent types of PLs (PE, LPC, . . ., see below). In this 
equation we have assumed that across each molecule's 
headgroup area the change in the local membrane cur- 
vature is so small, that the integration in ([T]) simplifies 
to a multiplication with the pivotal area, Apj, occupied 
by a lipid of type i. The pivotal area of lipid, Ap^i, is 
defined as the area that remains unchanged in its size 
upon spherical bending [2^. (See Fig. [3] for an illustra- 
tion.) Note that in writing ([2]) we disregard contributions 
due to orientational ordering [29j , and implicitly assume 
properly aligned lipids. 

We may construct the PL monolayer's average free 
Gibbs energy, G, per lipid as the sum of independent, sin- 
gle molecule energies, Ui, plus their corresponding con- 
figuration entropies, i.e. 

g G/N = - / dA^ [x,C/,(r) + x^fceTlna;,] , (3) 

with, Xi, the local lipid mole fraction in the membrane; 
n, the number of lipid species in the monolayer; and, N, 
the total number of lipids. In the special case of a planar 
membrane (H — > 0, AT — » 0) 5 approaches 

n 

9p ~ ^ ^ Xp^ifip^i, fip^i = 2fcc^iCQ j + k^T InXp^i, (4) 
1=1 

where we have used Xp^t and fip^i to denote the mole 
fraction and chemical potential of a flat membrane, re- 
spectively. 

According to ([3]), the elastic energy of a mixed mem- 
brane depends on both its composition and the local cur- 
vature of the surface. We are interested in whether or not 
the total, interfacial energy of curved membrane surfaces 
can be reduced by changing the lipid composition relative 
to the composition of the flat membrane. 

In equilibrium, the local lipid composition, Xi, in 
eq. ([3]) is given by the Boltzmann factor, 

exph([/,(r)-A^p,,)/fcBr] 
E;=iexph(C/,(r)-^p,,)/fcBT]- ^> 

and can be interpreted as the probability for a lipid of 
type i to populate regions of energy Ui{r), while the pla- 
nar part of the membrane acts as a reservoir, which ei- 
ther provides or absorbs lipids depending on the differ- 
ence between Xi and Xp^t. Note that for a flat membrane 
^ recovers the limiting planar lipid composition, i.e. 
{H -^Q,K ^Q)^x,^Xp,^. 
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FIG. 2: A LD is modeled as spherical cap of radius, _Rc, and 
height He, in an otherwise planar membrane (sectional view). 
Thick lines represent the cytosolic and luminal leaflets of the 
ER membrane. The two panels illustrate situations for differ- 
ent values of _Rc, but constant cap volume, Kap- 



A. Shape approximation 

Eq. ([3]) and ([5]) allow to estimate the energy change 
caused by local adaptations of both the lipid composi- 
tion and the surface curvature. Next, we make further 
assumptions regarding the shape evolution of emerging 
LD that form from initially flat ER membranes. 

In every stage of their biogenesis, nascent LD are as- 
sumed to form perfectly spherical protuberances of cap 
radius, i?c- The spherical cap approximation supposes a 
sudden transition from the planar to the spherical region 
of the membrane. 

Fig. [2] illustrates two possible geometrical configura- 
tions. In both cases the total volume of the LD, Vcap, is 
identical. In the following, it is convenient to discuss the 
effect of the LD volume in terms of the radius, i?LD, of 
an associated, full sphere of equal volume. 

Let i?c denote the radius of a spherical cap, then 
the cap surface Ac(i?c,J^LD) enclosing a cap volume 
14ap = 47ri?LD'^/3; as well as its corresponding pinch-off 



area Ai(i?c, i?LD), can be written as 



(6) 



Ac(i?c,i?LD) = 27ri?2i/(i?LD/i?c) 
^l(i?c,i?LD) = 7ri?2 [2i/(i?LD/i?c) - i/2(i?LD/i?c)], (7) 

'1 
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where i/c = RcH{Rlt:)/ Rc) gives the total height of the 
spherical calotte (see Fig. [2]). While i?LD is a measure of 
the LD volume, the radius Rc is the representative of the 
shape of the nascent LD. A LD detaches if its pinch-off 
area, Ai, vanishes. Thus for matured LD i?LD = ^^c and 
H{1) = 2, such that Ai ^ and A^ = A-kR^. 



B. Free Gibbs energy of a LD 

In the spherical cap approximation, the local and 
Gaussian curvature is independent of the position on 
the cap surface and simply given by H = 2/Rc and 
K = 1/R^, respectively. Thus, the free Gibbs energy 



5'cap(^c) = ^ [xiUi{Rc) + XikBTln: 



(9) 



depends parametrically on the cap radius. 

With an expression for the surface at hand it is possi- 
ble to estimate the average number of lipids on nascent 
LD as iVcap = X]r=i ^i^c/Ap^i, as well as the number of 
lipids, which originally were found in the pinch-off area 
but migrated as A^p = X]r=i ^p,i^i/^p.4- Then, the total 
Gibbs energy of formation for a LD becomes 

AG = AfcapScap - iVp5p. (10) 

In the following we will keep the total LD volume, Vcap, 
constant and ask which configuration, out of all possible 
values of Rc, is energetically favored. 



C. Parameter estimation 

Glycero-PL, in particular phosphatidic acid (PA), 
phosphatidyl choline (PC), PE, phosphatidyl inositol 
(PI), LPC and phosphatidyl serine (PS), are major con- 
stituents of ER membranes [2^. According to ((T|), our 
modeling approach requires four characteristic numbers 
per lipid species i. Two of which provide structural infor- 
mation {Ap^i and co,i), while the remaining two describe 
elastic properties (feci and fcc,i)- The latter were found 
to vary little for different lipid species [30| . and experi- 
mental results suggest that 
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irrespective of specific lipids [30| ■ 

However, literature values on geometrical data for 
lipids scatter substantially. Partly because data, such 
as the head-group area or chain-length, depend on tem- 
perature, salt content, the phase of the membrane, etc., 
but also because of variations in experimental methods. 
For instance, [sil ] presents ten published head-group ar- 
eas for fluid phase Dipalmitoyl-PC, ranging from 57 
to 73 A^. We have compiled representative values of 
structural data of lipids (Table |T| , being fully aware that 
these values may only be considered as rough estimates. 
A thorough discussion of our estimations may be found 
in the supplementary material lAl 

Here we describe lipids based on a geometrical pack- 
ing parameter, Si, which grossly characterizes the shapes 
of lipids by relating their entire volume Vi, to the vol- 
ume given by their head-group areas aoA, times the lipid 
length li, i.e. Si = Vi/ (ao.ili)- Then the spontaneous 
curvature, cg.i, of a monolayer of identical lipids for a 
cylindrical system can be expressed as [H, [s^ 
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In the last part of the equation we have approximated 
the volume of a PL by the volume of a truncated cone, 
Vi = (oo.i + bos)li/2, where ao,i and bo,i are the lipid's 
head- and base area, respectively (Fig. [3]). i?w.j, is typ- 
ically measured in fully hydrated lipid phases and refers 
to the radius of curvature of the lipid-water interface. 
For simplicity and despite substantial variations in the 
lipid's length, we assume that all lipid -typ es arc basi- 
cally as long as PE, i.e. Z = 22 A To further 
ease analysis we define the pivotal plane to sit midway 
between the base and head area of a lipid, thus 



R 



1/2. 



(13) 



However, according to experiments the neutral plane sits 
closer to the boundary between the hydrocarbons and 
polar group. For PE a representative value is 0.371 fs^ . 

A classification according to the packing parameter, 
S, reveals three major groups; PE and PA have 5*^ » 1, 
forming wedge-shaped, inverted truncated cones, while 
PC as well as PS, with S ~ 1, are quite cylindrical in 
shape. Finally, PI and LPC have S values between 0.7 
and 0.9, corresponding to truncated cones with lipid foot 
areas much smaller than their head-group areas. 

Finally, Table |n] lists PL composition of membranes 
for both, the LD and the ER [H. The two data 

sets agree well for PC and PE, showing a considerable 
reduction of both lipids in the LD membrane relative to 
the ER. Also, both data sets reflect a substantial increase 
in PA and PI. Conflicting results arc obtained for PS; 



TABLE I: Structural data of various PL. Literature values are 
marked by footnotes, estimated values are highlighted in gray, 
calculated values [using - (|13|l ] are printed without any 
tag, and framed lines contain data which enter the calculation. 
Abbreviations: ao, lipid head-group area; ^p, molecular area 
of the pivotal plane; bo, base area; 5", shape factor; co, spon- 
taneous curvature of the membrane monolayer; R„, cylinder 
radius to the lipid-solvent interface; Rp cylinder radius to the 
pivotal plane; fcc and kc, mean and Gaussian curvature elastic 
modulus, respectively. 
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FIG. 3: Characterization of different PL according to their 
structure factor, Si — {1 + bo,i/ao,i)/2, panel a and c. S > 1 
indicates cone like shape, which favors concave membranes. 
5 < 1 describes lipids as inverted cones, which produce convex 
curvatures. Panel b illustrates the pivotal plane, which by 
definition does not change in size upon spherical bending. 
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FIG. 4: Literature values for experimentally measured PL 
composition of ER and LD membranes compared to our ther- 
modynamic model at various radii R^. Data listed in 0, 
for "other lipids" have been redistributed such that each col- 
umn sums up to one. (For actual numerical values see sup- 
plemental material, Table |lTl) Note, that both outermost 
columns are identical. They have been duplicated to facili- 
tate easier comparison. 



while the fraction of PS is doubled in [23|, it is found to 
be slightly reduced in 0. Precise values for LPC are not 
known, but it has been reported that LPC is enriched in 
the LD membrane [i^ . To study the packing effects even 
on LPC, wc assume its x„,i equal to that of PA. 



5 



c = 1/Rc/[nm''] 

10 1 0.1 0.01 0.001 

■| L',AikibiLi^i..L.» ' ' ' ' ' ' ' ' - 




FIG. 5: Predicted lipid composition, Xi{Rc), of a PL mono- 
layer covering a spherical surface of radius Rc (lines in upper 
panel). For Rc oo, the curves slowly converge to the mea- 
sured composition Xp^i of the planar ER membrane [2^. The 
lipid composition of ripe, detached LD, as determined by [23 |. 
are indicated by symbols at the position of best fit (arrow). 
Lower panel, average free Gibbs energy per lipid, g as func- 
tion of cap radius, Rc- It shows the impact of spherical bend- 
ing on g for a compositionally optimized lipid membrane (full 
line) and a compositionally rigid membrane with frozen-in ER 
composition (dashed line). 



III. RESULTS 

The optimal PL composition, Xi, of a LD membrane 
was evaluated as function of cap radius, Rc, using 
Fig. [5] illustrates the predicted lipid distribution in the 
LD monolayer (upper panel) originatin g fr om a planar 
ER with PL composition according to [2J|. Rc repre- 
sents the radius of a spherical protuberance in an oth- 
erwise planar lipid layer (Fig. [2]). Its inverse, c = l/Rc, 
may be interpreted as the mean curvature of the associ- 
ated sphere and Fig. [5] can be understood in terms of the 
curvature dependence of the lipid composition. 

Fig. O describes the demixing process that occurs 
for curved membranes. For small curvatures (c < 
0.005 nm^^ ^ Rc > 200 nm), the local PL composi- 
tion in the ER membrane is only weakly affected and 
close to the lipid distribution of the unperturbed planar 
ER. On the other hand, for small cap radii {Rc < 100 nm 
<J4> c > 0.01 nm""'^ ), a significant demixing is observable; 
while generally lipids with Si > 1 (PE, PA, and PC; thin 
lines in Fig. [5]) decrease with decreasing Rc, the fractions 
of PS, PI, and LPC {Si < 1; thick lines) considerable in- 



crease. The observation is consistent with expectation, 
as lipids with S > I tend to form convex surfaces, while 
S < 1 favors concave topologies. However, this trend 
changes at high curvatures, where PI and PS content 
peak. Thus, for even higher curvatures, i.e. Rc < 1 nm, 
the monolayer enriches only on LPC at the expense of all 
other PL. 

PC is an exception to the rule stated above. Despite 
its S value being larger than 1 {Si = 1.11) it does not 
migrate but slightly increases for intermediate values of 
Rc- Eventually, however, even PC drifts off for i?c < 5 
nm. Only a small fraction of this discrepancy can be 
attributed to the differences in elastic moduli. For in- 
stance, if fcc,i is changed from 9kBT to lOfceT (plus the 
corresponding changes in kc^i), then the maximum PC 
content in the membrane decreases by only 2% (data not 
shown) . 

In the lower panel of Fig. [5] we show the change of 
Gibbs energy per lipid accompanied by the demixing 
process (full line) according to For comparison we 
have also plotted the Rc dependence for a membrane 
with frozen-in ER composition [dotted line; ^ together 
with the setting Xi = Xp^i,'iRc\. Both lines converge 
for Rc ^ oo. However, while for a membrane with 
fixed composition the energy increases with decreasing 
Rc, the behavior is more subtle for a membrane, which 
optimizes its lipid composition. In this case the Gibbs 
energy initially decreases with decreasing Rc, reaching a 
pronounced minimum at i?c ~ 9 nm. Thus an optimized 
packing is able to overcompensate the increase in elastic 
energy caused by spherical distortion of an originally pla- 
nar surface. Only at very high curvatures is regrouping 
of lipids unable to balance the rise in bending energy and 
the Gibbs energy sharply increases. 

To predict the total energy of formation for a nascent 
LD of volume AnRf^jy/S and cap radius Rc we evaluate 
PU)) . Fig.IHshows AG{Rc;Rld) as function of Rc. Note 
that AG depends parametrically on the LD-volume, rep- 
resented by i?LD- That is, each line in Fig.[B]shows the de- 
pendency of the energy on Rc for a fixed LD volume. De- 
pending on the size of the LD, that is on the value of Rhu, 
the free Gibbs energy shows either a single minimum at 
values Rc/Rld > 1 (for Rld < 4.2 nm), two minima - 
one at Rc/Rlb ~ 1 and the other at Rc/RhD > 1 (for 
4.2 nm < i?LD < 6.5 nm) -, or exactly one minimum at 
Rc/Rld = 1 (for Rld > 6.5 nm). 

Fig. [6] may be interpreted in such a way that for con- 
stant LD volume, minima in the Gibbs energy single 
out preferred configurations, i.e. identify likely values 
of Rc for a given LD volume (thick full and dashed line). 
Rc/Rld = 1 indicates that a LD docs no longer form a 
spherical cap, but a full sphere. Thus, if the minimum 
in energy shows up at such a point, the LD has been 
fully fabricated and detaches from the ER. If, however, 
a minimum appears at values Rc/Rld > 1, then the 
nascent LD forms a stable spherical protuberance of vol- 
ume 47ri?Lj-,/3 and cap radius Rc, but remains an integral 
part of the ER membrane. The position of the minima. 



6 




FIG. 6; Total free Gibbs energy, AG{Rc; Rld), as function 
of the normalized cap radius, Rc/Rw- Each full, thin line 
represents the energy dependence for a fixed but constant 
LD volume (TIld ~ [0.5, 10] nm, increment 0.5 nm). The 
observed (local) minima are marked by thick lines. For a 
better understanding of the meaning of the normalized x- 
axis the effective geometrical conformation is illustrated for 
two exemplary values on top. Rc denotes the cap radius, 
while i?LD represents the equivalent spherical radius of the 
cap. Inset: Position of the local energy minima, -Rc.opt, as 
function of LD size, _Rld. The thick lines correspond to the 
thick lines in the main figure. The dotted line indicate points 
where i?c,opt = Rho- 

-Rc,opt, as function of LD size, i?LD, is illustrated in the 
inset of Fig. [H] In this plot, points along the dotted line 
{Rc/Rld = 1) indicate fully detached LD. 

The stability of a nascent LD retained in the ER may 
be best estimated by the height of the energy barrier 
between the position of the minimum, i?c,opt, and the 
point where the cap becomes a full sphere, i.e. at Rc = 
RhD- This barrier shrinks with increasing LD volume. In 
our calculation we find that for LD diameters larger than 
2i?LD = 13 nm this barrier vanish completely and does 
not reappear again. Thus such LD are fully matured and 
detach from the ER. It is important to note that if we 
keep the lipid composition of the ER membrane fixed, 
then no decreasing energy barrier is observable (data not 
shown), and LD remain as integral part of the ER - in 
contrast to experimental findings. Thus lipid demixing 
is an essential part in the budding process. Moreover, we 
conclude that our lipid demixing model can explain the 
budding process based on biophysical reasons only. 

To understand the reason for the disappearance of the 
energy barrier consider the change in energy associated 
with LD formation, pO)) . In essence AG(i?c;-RLD) de- 
pends on the product between gcap(-Rc) and Ac{Rc', Rlb)- 
Since Ac{Rc; Rld) has values just for Rc > Rld, only 
these distances contribute to the product. On the other 



hand, gcap{Rc) does not depend on Rld, and neither does 
its minimum at i?c ~ 9 nm. Therefore, if Rld is small 
enough so that Rld lies to the left of the minimum in 
Fig. O then also the total energy exhibits a pronounced 
minimum and AG{Rc', Rld) rises sharply in-between. An 
energy barrier is formed. On the other hand, if Rld lies 
to the right of the minimum, then gcap(-Rc) is cut-off be- 
fore it reaches its minimum, a stabilizing energy barrier 
cannot form and LD detach from the ER. Thus the po- 
sition of the minimum in gca.p{Rc) in relation to the LD 
size, i.e. Rld , is decisive for the appearance of an energy 
barrier due to gc&piRc)- 

The predicted bud-off diameter depends on the initial 
ER membrane composition. To what extent the bud-off 
radius can be controlled by the ER membrane composi- 
tion has not been investigated. However, we successfully 
checked whether we get similar results if the PL compo- 
sition of LD rather than the one of the ER is used as an 
input for the simulation. Therefore we solved ([5]) such 
that at a typical LD size of Rc = 200 nm [1] the mea- 
sured PL composition [2^ is retrieved. In that case we 
predicted a bud-off diameter of 6.7 nm (see supplemen- 
tary materiaim Fig. [7] and Fig.!]). 

Using the parameter settings in Table [H our model 
is able to describe budding. At the bud-off diameter, 
2i?LD = 13 nm, we observed a significantly altered mem- 
brane composition in comparison with the ER. Consis- 
tently, we got a similar scale for the budding size, if we 
fit Rc such that the computed composition (taken from 
Fig. [5]) best approximates the measured data. By doing 
so, we obtained Rld = 18.5 nm and Rld = 3.2 nm for 
the data in [2^ and 0], respectively (see Fig. [3]). How- 
ever, in any case our predicted bud-off radius is an order 
of magnitude lower than the typical LD radius of 200 
nm Q . This suggested that if the standard model of LD 
biogenesis is correct, then it is a two-step process. Step 
one, growth of a nano-sized LD nucleus within the ER 
membrane. Step 2, cytosolic ripening of a nano-LD after 
they have budded-off from the ER. 



IV. DISCUSSION 

Lipid demixing in curved membrane regions. We have 
presented a thermodynamic model to describe LD forma- 
tion and investigated the interplay between membrane 
curvature and demixing in the PL composition. In fact, 
lipids favoring a convex surface {S < 1, i.e. PS, PI, 
LPC) are enriched while PE, PA, and PC are depleted. 
Curvature induced demi xing is hardly surprising and al- 
ready well documented [25|. All the more important, 
however, is the observation that the free energy per lipid 
for a curvature-optimized packing can drop below that 
of a planar monolayer. In other words, moderately de- 
forming an originally planar mixed PL monolayer does 
not necessarily cost energy, but can result in an energy 
gain - provided that the lipid system can optimize its 
composition. 
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Demixing forces induce membrane curvature. Our 
analysis showed that during hpid demixing of membranes 
energy is gained rather than consumed. Thus we argue 
that demixing at least supports if not drives membrane 
deformation. Such a view is evoked by Fig. [G] By increas- 
ing the NL content in the ER one moves between two 
isolincs (transition from 1 to 2'; note that due to normal- 
ization the transition appears as a tilted line even though 
i?c remained constant). At 2' the membrane perceives a 
"demixing force" which further distorts the monolayer 
and pushes the system toward a new minimum at point 
2, thereby sucking NL in the nascent LD. Note that in 
Fig. [6] demixing forces occur to the right of the thick line, 
while to its left forces are fueled by elastic energy. Thus 
demixing is a general feature of our model. It does not re- 
quire LD biogenesis proceeding along the energy minima 
in Fig. [6l 

At this point it has become clear that only the exis- 
tence of a spherical surface - as a topological alterna- 
tive to the plane - renders the molecular shape an en- 
ergetically relevant quantity. This makes it energetically 
worthwhile for the membrane to separate the different 
lipids according to their molecular shapes. In that sense, 
the spherical topology is the real agent of the predicted 
demixing process. 

Demixing - a •potential method to recruit proteins from 
the ER membrane? We speculate that not only lipids 
but also tracer components of the ER membrane are en- 
riched or depleted on nascent LD as dictated by their spe- 
cific molecular shape factor. In principle, such a mech- 
anism would be able to explain the clustering of spe- 
cific proteins on raised membrane domains. For example, 
cone-shaped membrane proteins (as long as they are not 
transmembrane proteins), should react to the positive 
curvature of developing LD. This mechanism is poten- 
tially able to explain the so far unaccounted accumula- 
tion of over-expressed caveolin-1 on LD [4l|,|43, check cer- 
melli2006]. Cavcolin-1 has a hairpin-like structure, with 
both terminal facing the cytosol. Thus, its structure fac- 
tor is smaller than 1 and therefore, like PS, enriches on 
the LD surface. Likewise, proteins with molecular shapes 
akin to PE should be depleted. This offers the interest- 
ing perspective of an efRcient, self-organized, biophysical 
"protein sorting" mechanism that performs the recruit- 
ing of proteins from the pool of ER proteins according to 
their molecular shape. 

The birth of a lipid droplet. We argued that PL 
demixing in developing LD is not a consequence of the 
membrane curvature, but rather its driving force. The 
existence of pronounced minima in the Gibbs energy 
(Fig.[6|) leads to demixing forces during LD growth. Thus 
the following thermodynamic, budding mechanism is sug- 
gested: A LD forms by accumulating NL within the 
ER membrane, thereby producing a membrane protu- 
berance. Its curvature is adjusted by a demixing force 
in such a way that it minimizes the interfacial energy of 
the cytoplasmic monolayer (see Fig. [HI transition 2' to 
2). As long as the NL volume, i.e. the size of the pro- 



tuberance measured by i?LD, is small enough, the LD is 
prevented from detaching by an insurmountable energy 
barrier. This barrier is traced back to the sharp rise in 
the elastic energy at very high curvatures. However, for 
NL volumes, V, for which i?LD = \/3V^/(47r) > 6.5 nm, 
the developing LD docs not buckle enough. The increase 
in bending energy is too small and a stabilizing energy 
barrier cannot mold. Thus, such large LD do not ex- 
ist within the ER membrane, as they have budded off 
already at an earlier stage of their development. 

At first it might appear counterintuitive that the elas- 
tic properties of a thin skin such as a PL monolayer 
should ultimately be responsible for the formation of LD. 
However, considering a LD suspension as some sort of 
macroemulsion, these ideas are far from new: in fact, it 
has been shown that type and stability of oil-in-water 
emulsions are determined by sign and value of the mono- 
layer's spontaneous curvature [43j . 

LD size distributions based on electron microscopic 
pictures for various yeast mutant strains, are well repro- 
ducible and rather narrow, showing values between 250 
nm and 550 nm in diameter, with a maximum at about 
400 nm Q . Nevertheless these radii are an order of mag- 
nitude larger than those predicted by our model. What 
could be the reason for this discrepancy? Obviously, the 
experimentally observed distribution reflects the size of 
the LD suspension in the cytosol. This distribution does 
not necessarily have to be identical to the size of freshly 
released LD; implying that there exists a second, subse- 
quent phase of LD ripening within the cytosol. Thus, 
small LD produced from the ER would undergo a series 
of fusions through conglomeration and coalescence until 
they reach their full-grown size Q. A number of obser- 
vations support this view. 

Although current experimental methods allow for the 
detection of 400 nm structures, a fully matured LD, 
which is still part of the ER, has never been observed 
in any cell. Neither have smaller ones. This is proba- 
bly because such LD are too small to be resolved with 
conventional experimental techniques |3l • Even with 
the usage of electron microscopy, which provides the ul- 
timate resolution, nascent LD would be difficult to ob- 
serve. By cutting a section of 80-90 nm through a fixed 
cell it is unlikely to section also an attached LD. Nev- 
ertheless, very recently large NL globules were detected 
by electron microscopy within special sub-compartments 
of the ER membrane in human hepatoma cells [45j . The 
authors hypothesized that tight binding of anomaly lipi- 
dated apolipoproteins suppressed detaching and allowed 
detecting giant NL structures within the ER membrane. 
Whether this observation presents the missing link in val- 
idating the standard model of LD formations remains to 
be seen. 

TAG accumulates either in LD dispersed in the cy- 
tosol or in membrane micro-domains [46|. These micro- 
domains are tiny TAG depots residing between two mem- 
brane leaflets [43] . Such depots were estimated to have a 
size between 22-28 nm which is roughly consistent with 
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our diameters. This would suggest that these small in- 
clusions are simply nascent LD. However, based on their 
data [4^, an ultimate conclusion is not possible. 

Artificially prepared PL membranes on a solid sur- 
face showed time-dependent changes in the membrane to- 
pography after externally initiating NL production [4^ . 
Droplet-like structures were formed at the surface with 
an average diameter of about 50 nm. Again, this is 
comparable to the size regime suggested by our model. 
Then, in a second process stage, smaller LD coalesced 
into larger structures. Additionally, LD fusion has also 
been shown in vivo (49l |. 

If LD biogenesis proceeds in a two step process, i.e. 
budding of nano-LD, which later merge, then fusion will 
dramatically change the surface to volume ratio and we 
also expect alterations in the PL composition of the LD, 
too [501 . In fact, if, for simplicity, we assume that LD 
bud-off at a diameter of 40 nm then 1,000 nano-LD arc 
needed to produce a typically sized LD. The total sur- 
face of these nano-LD will be ten times larger than the 
surface of the full-blown LD. This excess of PL may be 
used as a reservoir, which would allow the growing LD to 
re-mix its PL composition in order to accommodate the 
decreasing curvature. However, this reservoir will not be 
accessible as a whole, as it is unlikely that all 1, 000 nano- 
LD merge simultaneously in a single big fusion. Thus 
lipid re-mixing due to fusing of LD ultimately depends 
on the process of coalescence. Nevertheless we expect 
it to be less efficient then the original lipid demixing in 
the ER membrane. In fact, in our model we obtained 
PL compositions for stable nano-LD that are similar to 
those of fully matured ones (sec Fig. 2]). Thus, the lack 
of accessible PL-depots during the ripening process could 
explain why nano-LD and matured LD have similar com- 
positions, despite their different sizes. 

However, any remixing would change the composition 
but leaves the total number of excess PL unaltered. On 
the other hand, many LD are found to be surrounded 
by or adherent to double layer structures [sol . (unpub- 
lished, H Wolinski)]. These membranes are neither part 
of the ER nor do they form continuous structures with 
LD. Their origin and function is unclear. Here we spec- 
ulate that these membranes are a waste dump sprouting 
from excess PL during LD fusion. 

It is also interesting to link our results to mammalian 
very low-density lipoproteins (VLDL), which are another 
kind of lipid particles. Their NL core are of the same or- 
der of magnitude (30 — 90 nm) as our predicted nano-LD 
[Slj . This observation would support the so far untested 
hypothesis of ^] who suggested that ready-made nano- 
LD fuse with apoproteins to form VLDL. 

Model limitations. Our model assumes that LD, 
though forming between the leaflets of the ER membrane, 
lead to a protuberance only of the cytoplasmic leaflet, 
while the endoplasmic leaflet remains unaffected. Such 
an asymmetry may be real and caused by certain mem- 
brane proteins, differences in the PL composition of the 
leaflets or different properties of the solvents adjacent to 



the leaflets. However, a rigid endoplasmic leaflet is not 
an essential assumption. Our model could be adopted to 
account for LD growth in both leaflets. 

An essential part of our approach is the spherical cap 
approximation, which has been proven helpful in the con- 
text of the Young equation, which relates the contact an- 
gle of a liquid droplet on a solid surface to the inter-facial 
energies of all three phases involved [53j . However, this 
approximation neglects contribution from the neck re- 
gion. The neck connects the planar ER membrane with 
the LD monolayer and has the geometrical shape of a 
torus. In contrast to the cap, the neck region consist of 
both, positive and negative curvatures. 

Our current model is not able to address neck forma- 
tion. However, it was demonstrated theoretically that 
the neck shape is generated by an interplay between the 
local curvature and the membrane composition [H^, \EM ■ 
Similar to the observation reported here (accumulation 
of 5 < 1-lipids at high curvatures), it was shown that 
the neck region predominately consists of lipids with 
anisotropic, saddle-like intrinsic curvature. This is con- 
sistent with our results as these lipids are best adapted 
to the torical shaped neck. We are currently working 
on including the neck in our description. However, we 
do not expect any significant impact on our conclusions, 
because firstly the neck region is small in comparison to 
the spherical cap and secondly contributions to the total 
energy will be attenuated by lipid demixing. 



V. SUMMARY 

A biophysical model to describe LD biogenesis has 
been presented. We have estimated the local compo- 
sition of a PL monolayer, which is part of a spherical 
protuberance in an otherwise planar membrane. Four 
major results are obtained, (i) In order to use the avail- 
able surface area economically, PL with shape factors, 
5 < 1 (PI, PS, LPC), tend to accumulate on spherical 
surfaces. Conversely, lipids with S > 1 (PE, PC, PA) 
migrate. Therefore, a spherical protuberance of radius, 
_Rc, in an otherwise planar membrane causes a local lipid 
demixing on the curved membrane monolayer, (ii) As a 
consequence of the optimized packing of lipids, the energy 
per lipid gcap(^c) drops with decreasing Rc, i.e., with in- 
creasing curvature. Thus, induced by lipid demixing, a 
monolayer gains energy by forming a spherical protuber- 
ance of weak to intermediate curvature. Only at very 
high curvatures, gcapiRc) increases again, leading to a 
clear minimum in gcap(-Rc) at 9 nm. (iii) This energy 
minimum is key in understanding the formation of an 
energy barrier, which controls the budding of LD. (iv) 
The height of the energy barrier depends on the volume 
of the protuberance; it stabilizes protuberances of smaller 
volumes, but completely vanishes at higher volumes. Our 
model predicts that LD detach from the ER at a diameter 
in the order of 13 nm. 

In this paper we have suggested that LD formation 
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is driven by lipid dcmixing. This finding is based on 
a model calculation. The main ingredient entering this 
model is the geometrical structure of various lipids, de- 
scribed by the shape factor. Thus, by changing the 
shape factor, we expect different demixing effects and 
ultimately different bud-off sizes, which could be verified 
by m vitro measurements, similar to those performed by 
[48|. Experimentally the shape factor may be influenced 
by varying the chain length of FA in lipids, or by manip- 
ulating the pH of the solvent which influences the head 



group areas. 
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APPENDIX A: ESTIMATING STRUCTURAL 
PARAMETERS OF PL 

PE is the only lipid for which both, ao,j; = 54 and 
i = 105 have been determined experimentally [s^ . 
Together with cq,, = -0.043 A'^ [H] and dH]), one ob- 
tains I = 22 A, which is roughly half the thickness of a 
typical PE bilayer ranging between 40 and 50 A [lil. It 
also compares well with the value / = 18 A obtained by 

To estimate data for PC, we followed an idea 
of Marsh [sO], who estimated cg values for lipid 
dioleoylphosphatidyl-PE/DOPC mixtures by fitting ex- 
perimental data [28[. Wc assumed oo.i of PC to be 72 
A^ fsil . [3^ and used (fT2|) as fitting function for the data 
in (28[. An excellent fit is obtained, giving a value of 
bo,i = 87.8 A^ and subsequently cn ;, — -0.01 A~\ which 
agrees with the values given in [sl, HI, |3§| . 

For PA, PI, and PS the situation is less clear. As we 
have been unable to find experimentally measured head- 
group areas of PA and PI, we have estimated these values 
by extrapolating the linear relationship between molec- 
ular volume and head-group area obtained for PC, PE 
and 0-LPC. Also, the molecular area at the pivotal plane 
Ap^i are unknown for PS, PI and PA; again, values have 
been guessed by comparison with corresponding values 
of other lipids. 



TABLE II: Literature values for experimentally (E) measured 
PL composition of ER (xp^i) and LD (xi) membranes com- 
pared to the results of our thermodynamic model (C). Data 
listed in 0, for "other lipids" have been redistributed such 
that each column sums up to one. Rc values denote the radii 
of best fit for the experimental data. Calculated PL compo- 
sitions represent values at the radii of best fit. 









x^{C) 


x,{E') 


x^{C) 


PE 


0.336 


0.243 


0.266 


0.203 


0.059 


PC 


0.515 


0.382 


0.549 


0.368 


0.537 


PS 


0.068 


0.141 


0.083 


0.055 


0.192 


PI 


0.077 


0.208 


0.097 


0.32 


0.202 


O-LPC 


0.002 


0.013 


0.003 


0.027 


0.009 


PA 


0.002 


0.013 


0.002 


0.027 


0.001 


-Rc (nm) 


00 


18.5 


18.5 


3.2 


3.2 



'Ref. m ''Ref. [1 
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APPENDIX B: CHANGING THE PL 
COMPOSITION OF THE ER MEMBRANE 



c = 1/Rc/[nm"'] 
0.1 




1000 



FIG. 7: Predicted lipid composition, Xi{Rc), of a PL mono- 
layer covering a spherical surface of radius Rc (lines in upper 
panel). The PL composition of the ER membrane was ad- 
justed such that at _Rc = 200 nm (arrow) the measured PL 
composition of a LD is met [2^. 




FIG. 8: Total free Gibbs energy, AG{Rc; -Rld), - for the data 
of Fig. [7]- as function of the normalized cap radius, Rc/Rld- 
Each full, thin line represents the energy dependence for a 
fixed but constant LD volume (i?LD = [0.5, 10] nm, incre- 
ment 0.5 nm). The observed (local) minima are marked by 
thick lines. For a better understanding of the meaning of the 
normalized s-axis the effective geometrical conformation is il- 
lustrated for two exemplary values on top. Inset: Position of 
the local energy minima, i?c,opt, as function of LD size, _Rld. 
The thick lines correspond to the thick lines in the main fig- 
ure. The dotted line indicate points where -Rc,opt = Rlu- 



TABLE III: Experimentally (E) and computationally (C) ob- 
tained PL composition for ER (xp,i) and LD (xi) membranes. 
The calculated PL composition of the ER membrane was ad- 
justed such that at Rc = 200 nm the measured PL composi- 
tion of a LD is met [S^. Data listed in [2^ for "other lipids" 
have been redistributed such that each column sums up to 



one. SrX = 




— x]/x 


denotes the relative 


error with 


respect to x. 














x,{E) 


Xp,i{E) 


^P.i(C) 


SrXp^i 


SrXi 


PE 


0.243 


0.336 


0.248 


-0.262 


0.021 


PC 


0.382 


0.515 


0.381 


-0.260 


-0.003 


PS 


0.141 


0.068 


0.140 


1.059 


-0.007 


PI 


0.208 


0.077 


0.205 


1.662 


-0.014 


0-LPC 


0.013 


0.002 


0.013 


5.5 


0.000 


PA 


0.013 


0.002 


0.013 


5.5 


0.000 


Rc (nm) 200 


oo 


oo 







